home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / slaqge.z / slaqge
Text File  |  1996-03-14  |  3KB  |  133 lines

  1.  
  2.  
  3.  
  4. SSSSLLLLAAAAQQQQGGGGEEEE((((3333FFFF))))                                                          SSSSLLLLAAAAQQQQGGGGEEEE((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SLAQGE - equilibrate a general M by N matrix A using the row and scaling
  10.      factors in the vectors R and C
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE SLAQGE( M, N, A, LDA, R, C, ROWCND, COLCND, AMAX, EQUED )
  14.  
  15.          CHARACTER      EQUED
  16.  
  17.          INTEGER        LDA, M, N
  18.  
  19.          REAL           AMAX, COLCND, ROWCND
  20.  
  21.          REAL           A( LDA, * ), C( * ), R( * )
  22.  
  23. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  24.      SLAQGE equilibrates a general M by N matrix A using the row and scaling
  25.      factors in the vectors R and C.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      M       (input) INTEGER
  30.              The number of rows of the matrix A.  M >= 0.
  31.  
  32.      N       (input) INTEGER
  33.              The number of columns of the matrix A.  N >= 0.
  34.  
  35.      A       (input/output) REAL array, dimension (LDA,N)
  36.              On entry, the M by N matrix A.  On exit, the equilibrated matrix.
  37.              See EQUED for the form of the equilibrated matrix.
  38.  
  39.      LDA     (input) INTEGER
  40.              The leading dimension of the array A.  LDA >= max(M,1).
  41.  
  42.      R       (input) REAL array, dimension (M)
  43.              The row scale factors for A.
  44.  
  45.      C       (input) REAL array, dimension (N)
  46.              The column scale factors for A.
  47.  
  48.      ROWCND  (input) REAL
  49.              Ratio of the smallest R(i) to the largest R(i).
  50.  
  51.      COLCND  (input) REAL
  52.              Ratio of the smallest C(i) to the largest C(i).
  53.  
  54.      AMAX    (input) REAL
  55.              Absolute value of largest matrix entry.
  56.  
  57.      EQUED   (output) CHARACTER*1
  58.              Specifies the form of equilibration that was done.  = 'N':  No
  59.              equilibration
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSLLLLAAAAQQQQGGGGEEEE((((3333FFFF))))                                                          SSSSLLLLAAAAQQQQGGGGEEEE((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              = 'R':  Row equilibration, i.e., A has been premultiplied by
  75.              diag(R).  = 'C':  Column equilibration, i.e., A has been
  76.              postmultiplied by diag(C).  = 'B':  Both row and column
  77.              equilibration, i.e., A has been replaced by diag(R) * A *
  78.              diag(C).
  79.  
  80. PPPPAAAARRRRAAAAMMMMEEEETTTTEEEERRRRSSSS
  81.      THRESH is a threshold value used to decide if row or column scaling
  82.      should be done based on the ratio of the row or column scaling factors.
  83.      If ROWCND < THRESH, row scaling is done, and if COLCND < THRESH, column
  84.      scaling is done.
  85.  
  86.      LARGE and SMALL are threshold values used to decide if row scaling should
  87.      be done based on the absolute size of the largest matrix element.  If
  88.      AMAX > LARGE or AMAX < SMALL, row scaling is done.
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.